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Abstract. The unavoidable irreversible loss of power in a heat engine is found 
to be of quantum origin. Following thermodynamic tradition a model quantum 
heat engine operating in an Otto cycle is analyzed, where the working medium 
is composed of an ensemble of harmonic oscillators and changes in volume 
correspond to changes in the curvature of the potential well. Equations of motion 
for quantum observables are derived for the complete cycle of operation. These 
observables are sufficient to determine the state of the system and with it all 
thermodynamical variables. Once the external controls are set the engine settles 
to a limit cycle. Conditions for optimal work, power, and entropy production are 
derived. At high temperatures and quasistatic operating conditions the efficiency 
at maximum power coincides with the endoreversible result ijq = 1 — \/Tc/Th- 
The optimal compression ratio varies from C = ^jTh/T c in the quasistatic limit 
where the irreversibility is dominated by heat conductance to C = (T^/Tc) 1 ^ 
in the sudden limit when the irreversibility is dominated by friction. When the 
engine deviates from adiabatic conditions the performance is subject to friction. 
The origin of this friction can be traced to the noncommutability of the kinetic 
and potential energy of the working medium. 
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1. Introduction 

Our cars, refrigerators, air- conditioners, lasers and power plants are all examples of 
heat engines. Practically all such heat engines operate far from the ideal maximum 
efficiency conditions set by Carnot 1 . To maximize the power output efficiency 
is sacrificed. This trade-off between efficiency and power is the focus of finite time 
thermodynamics" . The field was initiated by the seminal paper of Curzon and Ahlborn 
[2]. From everyday experience the irreversible phenomena that limits the optimal 
performance of engines [3] can be identified as losses due to friction, heat leaks, and 
heat transport. Is there a unifying fundamental explanations for these losses? Is it 
possible to trace the origin of these phenomena to quantum mechanics? To address 
these issues we follow the tradition of thermodynamics constructing a model quantum 
heat engine. 

Gedanken heat engines are an integral part of thermodynamical theory. Carnot 
in 1824 set the stage by analyzing an ideal engine Q]|I|- Carnot's analysis preceded 
the systematic formulation that led to the first and second laws of thermodynamics 

Amazingly, thermodynamics was able to keep its independent status despite the 
development of parallel theories dealing with the same subject matter. Quantum 
mechanics overlaps thermodynamics in that it describes the state of matter. But 
in addition, quantum mechanics includes a comprehensive description of dynamics. 
This suggests that quantum mechanics can originate a concrete interpretation of the 
word dynamics in thermodynamics leading to a fundamental basis for finite time 
thermodynamics . 

The following questions come to mind: 

• How do the laws of thermodynamics emerge from quantum mechanics? 

• What are the quantum origins of irreversible phenomena involving friction and 
heat transport? 

• What is the relation between the quasistatic thermodynamical process and the 
quantum adiabatic theorem? 

To address these issues a model of a reciprocating quantum heat engine is constructed. 
Extreme care has been taken to choose a model which can be analyzed from first 
principles. The specific heat engine chosen is a quantum version of the Otto cycle. 

The present study follows a tradition of analyzing quantum models of 
reciprocating heat engines El EI El HD1 HH H21 H31 H3| US]- Most of these models 
can be classified as endoreversible. In such models the source of irreversibility can 
be traced to the heat transfer between the working medium and the heat baths. 
Internal irreversibility was studied in a coupled spin working medium . The 
present paper employs an ensemble of harmonic oscillators as a working medium which 
incorporates both internal and external sources of irreversibility. A quantum version 
of the Otto cycle is described in section [21 The link between the quantum state and 
thermodynamical observables is developed in section [3J This model allows a closed 
form solution of the quantum dynamics (section 0J) . At low temperature the results 
have similarities with the spin model, at higher temperatures the well studied finite 
time thermodynamical performance characteristics appear naturally (sectional). 
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2. The quantum Otto cycle 

Nicolaus August Otto invented a reciprocating four stroke engine in 1861 and won 
a gold medal in the 1867 Paris world fair JHj- The basic components of the engine 
include hot and cold reservoirs, a working medium, and a mechanical output device. 
The cycle of the engine is defined by four branches: 

(i) The hot isochore: heat is transferred from the hot bath to the working medium 
without a volume change. 

(ii) The power adiabat: the working medium expands producing work while isolated 
from the hot and cold reservoirs. 

(hi) The cold isochore: heat is transferred from the working medium to the cold bath 
without a volume change. 

(iv) The compression adiabat: the working medium is compressed consuming power 
while isolated from the hot and cold reservoirs. 

Otto already determined that the efficiency n of the cycle is limited to r\ Q < 

1 — (^ L )~ _1 (where V c /h and T c / h are the volume and temperature of the working 
medium at the end of the hot and cold isochores respectively) 0. As expected, the 
Otto efficiency is always smaller than the efficiency of the Carnot cycle r\ < r\ c = 



2.1. Quantum dynamics of the working medium 

The quantum analogue of the Otto cycle requires a dynamical description of the 
working medium, the power output and the heat transport mechanism. 

The working medium is composed of an ensemble of non interacting particles 
described by the density operator p. These particles are confined by an harmonic 
potential V(Q) — \Q 2 ■ Expansion and compression of the working medium is carried 
out by externally controlling the curvature of the potential k(t). The energy of an 
individual particle is represented by the Hamiltonian operator: 

" - I/ 2 + tt« 2 • <» 

where m is the mass of the system and P and Q are the momentum and position 
operators. All thermodynamical quantities will be intensive i.e. normalized to the 
number of particles. In the macroscopic Otto engine the internal energy of the working 
medium during the adiabatic expansion is inversely proportional to the volume. In 
the harmonic oscillator, the energy is linear in the frequency uj(t) — ^k(t)/m |17| . 
This therefore plays the role of inverse volume . 

The Hamiltonian is the generator of the evolution on the adiabatic branches. 
The frequency w changes from LUh to lu c in a time period The in the power adiabat 
{ujh > w c ) and from ui c to u>h in a period T c h in the compression adiabat. The dynamics 
on the state p during the adiabatic branches is unitary and is the solution of the 
Liouville von Neumann equation |18|: 

P(t) = -~{H(t),p(t)\ > (2) 

where H is time dependent during the evolution. Notice that [H(t), H(t')] ^ since 
the kinetic energy does not commute with the varying potential energy. 
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The dynamics on the hot and cold isochores is an equilibration process of the 
working medium with a bath at temperature Th or T c . This is the dynamics of 
an open quantum system where the working medium is described explicitly and the 
influence of the bath implicitly |191 12U| : 

Ht) = - l -[H,p] + C D {p) , (3) 

where Cu is the dissipative term responsible for driving the working medium to 
thermal equilibrium, while the Hamiltonian H — H{ui) is static. The equilibration is 
not complete since only a finite time t>> or r c is allocated to the hot or cold isochores. 
The dissipative "superoperator" Cd must conform to Lindblad's form for a Markovian 
evolution ^] , and for the harmonic oscillator can be expressed as |21l [221 : 

£d(p) = ki{a)pa - ]-{aa),p}) + k^{apa) - -{a<a,p}) . (4) 

where {A, B} = AB + BA. fej and fcj, are heat conductance rates obeying detailed 

k — ,,4- ^ 

balance rr = e k b T , where T is either Th or T c . a T and a are the raising and 

K x 

lowering operators. Notice that they are different in the hot and cold isochores since 

3 = ^VWQ + *\[^L P ) Spends on co. 

Equation Q is known as a quantum Master equation. It is an example of a 
reduced description where the dynamics of the working medium is sought explicitly 
while the baths are described implicitly by two parameters: the heat conductivity 
r = k[ — fcf and the bath temperature T. The Lindblad form of (0J guarantees 
that the density operator of the extended system (system+bath) remains positive (i.e. 
physical) ^Hj- Specifically, for the harmonic oscillator, Q has been derived from first 
principles by many different authors |22| . 

To summarize, the quantum model Otto cycle is composed of a working fluid 
of harmonic oscillators l[T)). with volume corresponding to 1/lu. The power stroke is 
modelled by the Liouville von Neumann equation J2} while the heat transport via 
a Master equation (l.'il 11) . It differs from the thermodynamical model in that a finite 
time period is allocated to each of these branches. Solving these equations for different 
operating conditions allows to obtain the quantum thermodynamical observables. 

3. Quantum thermodynamics 

Thermodynamics is notorious in its ability to describe a process employing an 
extremely small number of variables. A minimal set of quantum expectations {X n } 
constitutes the analogue description where (X n ) — Tr{X n p}. The dynamics of this 
set is generated by the Heisenberg equations of motion 

X = ^ + l -[H,X]+C* D {X) (5) 

where the first term addresses an explicitly time dependent set of operators, X(t). 

To move from dynamics to t/iermodynamics we need to consider the three laws. 
The energy expectation E is obtained when X = H , i.e. E = (H). The quantum 
analogue of the first law of thermodynamics |231 124] dE = dW + dQ is obtained by 
inserting H into JSJ: 

E = W + Q = ( — ) + (C* D (H)) ■ (6) 
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The power is identified as V = W = ( ^M- ). The heat exchange rate becomes 
Q = ( C* D (H) ) . The Otto cycle contains the simplification that power is produced or 
consumed only on the adiabats and heat transfer takes place only on the isochores. 

The third law is trivially manifest for a system at thermal equilibrium (Gibbs 
state) at temperature zero. The intricacies of proving the second law in the general 
case are beyond this paper. Instead, in accordance with thermodynamic tradition, its 
applicability will be demonstrated for a closed cycle. Note that thermodynamics will 
be maintained even though the system is not generally in thermal equilibrium. 

The thermodynamic state of a system is fully determined by the thermodynamical 
variables. Statistical thermodynamics adds the prescription that the state is 
determined by the maximum entropy condition subject to the constraints set by the 
thermodynamical observables [251 1261 127] . Maximizing the Von Neumann entropy 18 

Svn = -k B Tr{p\n{p)} (7) 

subject to the energy constraint leads to thermal equilibrium |27j . 

Peg = \ er ^ ( 8 ) 

H 

where fee is Boltzmann constant and Z = Tr{e k B T } is the partition function. 

The state of the working medium in general is not in thermal equilibrium. In order 
to generalize the canonical form © additional observables are required to define the 
state of the system. The maximum entropy state subject to this set of observables 
(Xj) = tr{Xjf>} becomes 

P = |cxp \T > ( 9 ) 

where j3j are Lagrange multipliers. The generalized canonical form of is meaningful 
only if the state can be cast in the canonical form during the entire cycle of the engine 
leading to j3j — @j(t). This requirement is called canonical invariance |28|. A necessary 
condition for canonical invariance is that the set of operators X in ijjj} is closed under 
the dynamics generated by the equation of motion. If this conditions is also sufficient 
for canonical invariance, then the state of the system can be reconstructed from a small 
number of quantum observables (Xj)(t), which are the thermodynamical observables 
in the sense that they define the state under the maximum entropy principle. 

The condition for Canonical invariance on the unitary part of the evolution taking 
place on the adiabats is as follows. If the Hamiltonian is a linear combination of the 
operators in the set H(t) = ^2 m h m X m (h m (t) are expansion coefficients), and the 
set forms a closed Lie algebra [Xj,Xk] — J2i C\ k Xi (where C/ fe is the structure factor 
of the Lie algebra), then the set X is closed under the evolution [221 ■ I n addition 
canonical invariance prevails |3(J| . 

For the Otto cycle the set of the operators P 2 , Q 2 and D = (QP + PQ) form 
a closed Lie algebra. Since the Hamiltonian is a linear combination of the two first 
operators of the set (P 2 and Q 2 ), canonical invariance will result on the adiabatic 
branches. 

On the isochores the set has to be closed also to the operation of Cd- The set 
P 2 , Q 2 and D is also closed to Cd defined by l@J. This condition is only necessary 
for canonical invariance. Nevertheless for the harmonic working medium and Cd 
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defined in the closure is also sufficient for canonical invariance to take place (Cf. 
|Appendix A| ). 

The significance of canonical invariance is that all quantities become functions 
of a very limited set of thermodynamic quantum observables (Xj). The choice of 
operators Xj should reflect the most telling thermodynamical variables. Explicitly for 
the current engine the thermodynamical variables are chosen to be the expectations 
of the following operators: 

• The Hamiltonian H(t) = ^P 2 + \mu(t) 2 Q 2 . 

• The Lagrangian L(t) = ^P 2 - \mu(t) 2 Q 2 . 

• The position momentum correlation D = (QP + PQ). 

These operators are linear combinations of the same Lie algebra as Q 2 , P 2 and D. 

To make explicit the connection between the expectation values of the closed set 
of operators and the state of the system p an alternative product form for the density 
operator is used. It is defined by the parameters /3, 7 and 7*: 

p = le^V^e 7 *^ 2 , (10) 

where H = ^(aat + at a), D = -ih(a 2 - of) , L = -?f(a 2 + at*) and 

e 2 

From H1(J|I the expectations of H and a 2 are extracted leading to 

~\ hu(e 2 ^ - 4 77 * - 1) 

/ 2((e^ - l) 2 - 4 77 *) [ ' 

and 

2 7 * 



( e 0Ku _ !)2 _ 4 77 * • 
Equations l|12|) and l|13fl can be inverted leading to 

Huj(2(l) +iuj(b 



(13) 



7 = — 2 —^2 — ( 14 ) 

and 

e»» = < i y-<x) 2 +^ 2+ (°) 2 \ 2 (15) 

4^i\ +lu 2 (d^ - - 2 

Equations Ijl4(l and (|15(l make explicit the relation between the state of the system p 
(by equation ifTHjl ) and the thermodynamical observables (H) ,(D) and (L). 
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3.1. Entropy balance 

In thermodynamics the entropy S is a state variable. Shannon introduced entropy as 
a measure of the missing information required to define a probability distribution p 
|31| . The information entropy can be applied to a complete quantum measurement of 
an observable represented by the operator O with possible outcomes pj : 

S Q = ~kB^2 Pj 1ii Pj > ( 16 ) 

3 

where pj = Tr{Pjp}. The projections Pj are defined using the spectral decomposition 
theorem O = Y^j XjPj, where Xj are the eigenvalues of the operator O. S Q is then the 
measure of information gain obtained by the measurement. The Von Neumann entropy 
is equivalent to the minimum entropy S Y associated with a complete measurement of 

the state p, by the observable Y n where the set of operators Y n includes all possible non 
degenerate operators in the Hilbert space. The operator that minimizes the entropy 
commutes with the state [p, Y m i n ] = 0. Obviously Svn < Sq. This supplies the 
interpretation that Svn is the minimum information required to completely specify 
the state p. 

The primary thermodynamic variable for the heat engine is energy. The entropy 
associated with the measurement of energy Se = S^ in general differs from the Von 
Neumann entropy Se > Svn- Only when p is diagonal in the energy representation, 
such as in thermal equilibrium JSJl, Se — Svn- 

The Von Neumann entropy is invariant under a unitary evolution [30] . This is 
the result of the property of unitary transformations where the set of eigenvalues 
of p' = UpW is equal to the set of eigenvalues of p. Since the von Neumann 
entropy Svn is a functional of the eigenvalues of p it becomes invariant to any unitary 
transformation. In particular when applying a unitary transformation generated by 
the Hamiltonian H' to a thermal state p(0) = exp(- H/k b T)/Z then if [H f , p(Q)] ^ 0, 
the unitary transformation will always increase the energy entropy, > Se (0). 

Nevertheless the Von Neumann entropy will stay invariant Sv jv {t) = Sv n (0) = Se (0) . 

To find the Von Neumann entropy for the oscillator, the product-form of the 
density operator i|10|) is rewritten as 

Z 

where the relation between the coefficients \ in l|17fl and 7 and (3 in the product form 
IjlUfl are evaluated in |Appcndix B| , equation (|B.5I) . 

Employing (|17fl the Von Neumann entropy becomes: 



1 /^\ „ , , 2 



Svn = -X2 H +Im( X i)- D +Re ( Xl )— L + In - v . L Q A) 



csch 



( §hw 

V 2 



Equation l|18f) shows that Svn is a functional of the thermodynamical variables (H) , 
(D) and (L). 

The energy entropy Se of the oscillator (not in equilibrium) is found to be 
equivalent to the entropy of an oscillator in thermal equilibrium with the same energy 
expectation: 



Sb - W + t H L HM /-v -( 19 ) 
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Se in (|19fl is completely determined by the energy expectation E = \Hy. As an 

extreme example, for a squeezed state Svn — and Se > 0. 

In a macroscopic working medium, the internal temperature can be defined 
from the entropy and energy variables 1/Tmt = {§e~)v a * cons t an t volume. For 
the quantum Otto cycle, Se is used to define the inverse internal temperature 
I/Tint = (^e~) uj - Tint is a generalized temperature appropriate for non equilibrium 
density operators p. Using this definition the internal temperature Tj nt of the oscillator 
working medium can be calculated implicitly from the energy expectation: 

which is equivalent to the relation between temperature and energy in thermal 
equilibrium. 

The entropy expressions (|18|1 and Ijl9(l for the working medium illustrate that 
entropy changes along the cycle trajectory. Once the cycle reaches a limit cycle (Cf. 
Sec. (|5.1|l 'l then the internal entropy changes become cyclic. As a result any entropy 
change becomes zero |1U|. The irreversible character of the operation concentrates 
therefore in the entropy production in the thermal baths. Assuming a structureless 
reservoir, the entropy production is calculated from the heat transfer from the working 
medium: 

AS U = AS h + AS C = - ^- f " Q h dt - A. f ' Q cdt . (21) 

J-h JO ±c Jo 

where Th and T c are the hot and cold bath temperatures and Qh Qc are the heat 
fluxes Cf. equation JBJ. 



4. Solving the dynamics of the quantum Otto cycle 

The strategy to solve the dynamics of the working fluid along the engine's cycle of 
operation relies on obtaining closed form equations for the quantum thermodynamical 
observables and on using them to reconstruct the state of the system using the 
maximum entropy relation 



4-1- Heisenberg dynamics on the Isochores 

The working medium dynamics on the isochores represent an approach to thermal 
equilibrium where the Hamiltonian H is constant (u> — oj h / c is constant). The 
Heisenberg equations of motion for an operator X are 

X = =[H,X] +k l (a t Xa - -{at a, X}) + fcj {aXa) - ]^{dd\ X}) (22) 

Equation (|22H is the analogue of © and i n the Schrodinger frame. 

For the dynamical set of observables, the equations of motion become 







( 


d 


L 


(*) = 

V 


dt 


b 
\ i J 



-r o 

o -r 

o 4 

o o 



-u: 2 

-r o 
o o 
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\ i 


) 



(t) (23) 
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where V — k± — k-t is the heat conductance and fcf/fci = e - hu] / k BT 0De y S detailed 
balance where to = uj h / c and T = T h / C are defined for the hot or cold bath respectively. 
From 10 the heat current can be identified as: 

Q=-r m -(H }eq ) = r-(coth(^)-coth(^-)j24) 

where Tb is the bath temperature. In the high temperature limit the heat transport 
law becomes Newtonian: Q — T(Tb — Tint). 

The solution of the isochore dynamics (|23|) becomes: 

H(t) = e- rt (H(0) - I{H) eq ) + I{H) eq (25) 

( L \ Tt ( cos (2crf) -f sin (2ut) \f L\ 

{ D ) [ ) e { I sin (2wi) cos (2^) J ^ £> J W ^ 

where the equilibrium value of the energy is E = (H) eq = 4=r coth( 2 ^ T ) ■ 

The equilibration dynamics are characterized by the energy exponentially relaxing 
to its equilibrium value. The Lagrangian L and the momentum-position correlation 
D show a damped oscillation with a frequency of 2uj to an equilibrium value of zero. 
The identity operator I becomes a constant of motion representing the conservation 
of norm. 



4-2. The dynamics on the Adiabats 



On the adiabats the oscillator frequency changes from lu^ to lu c on the power expansion 
branch and from lo c to to^ on the compression branch. As a result the Hamiltonian 
is explicitly time dependent and [Hit), H (t')] ^ 0. Under these conditions the 
Heisenberg equations of motion JSJ for the dynamical set of operators become 



d_ 

dt 







( 


UJ 
UJ 


UJ 
UJ 


L 
D 


(*) = 




UJ 
UJ 




UJ 
uJ 

4 


V / J 















2 

-LI 






\ 






L 




D 


/ 


{ i J 



(t) 



(27) 



In general all operators in 1)27(1 are dynamically coupled. The coupling can be 
characterized by the nonadiabatic parameter a = It is important to notice that 
a has no explicit dependence on the mass of the oscillator. This means that the 
nonadiabatic character is universal. For example from the first law © the power can 
be identified as 

V = a ((£) - (X)) . (28) 

The power on the adiabats l|28|) can be decomposed to the "useful" external power 
Vex — Oi{H) and to the power invested to counter friction Vf = —a{L). It is "friction" 
in the sense that it reduces the power and is zero when the motion is infinitely slow 
and increases with speed (see below). The translation of this effect into the actual 
reduction of power, however, cannot be separated from the dissipative branches. 

A solution of the dynamics depends on an explicit dependence of ui(t) on time, 
and it can always be carried out numerically. An illustrative example is found when 
the non-adiabatic parameter a = — is constant during the adiabatic branches. This 

/ \t/r hc 

leads to the explicit time dependence ut(t) — ujqc = uj h I ^ J on the power 
adiabat, and a similar expression for the compression adiabat. Under these conditions 
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(29) 



the time derivative in (|27[) becomes stationary allowing a closed form solution to be 
obtained by diagonalizing the matrix in l|27|). 

Two simplifying limits exist. The first is the quasistatic limit, which is defined 
by a -C u>. In this case the fast motion (neglecting a) can be evaluated first in (|27ll 
leading to a zeroth order expression for L and D. From that expression one can 
construct an ansatz correct to first order in — : 

L(t) = ^i(0)co S (^)-o;(t)(^(0) + |^(0)) S in(^ 

= -^#(0) + (%H(0) + D(pj) cos(0) + £L(0) sin(0) 

where (pit) — J* 2uj(t)dt = 2 ("W~"°) j s the accumulated phase. For small ^, <f> diverges 
leading to fast oscillation of L[t). As a result L can be averaged out in the equation 
of the slow variable H leading to: 

^=*H (30) 

with the solution H(t) = ^^-H(ujo = uih) f° r the power adiabat and H(t) = 

^p-H(ujQ — ui c ) for the compression adiabat. 

The outcome of the quasistatic conditions is that the Hamiltonian follows 
adiabatically the change in frequency u>. This limit is equivalent to the invariance 
of the number operator N = j^H — \l in the evolution, i.e. (N) = constant. In this 
limit the power invested to counter friction is zero, Vf = 0. 

A higher order expression for the energy H can be obtained by integration, and 
will be useful in the consideration of friction. To second order in — the energy H on 
the adiabat becomes 

8uj%uj(t)H{t) = H(0) (-a 3 w sin(0) - 2a 2 Lu u(t) cos(</>) + 2a 2 u(t) 2 + 8u 2 u 

+ L(0) (-Aau 2 uj(t) sin(0) + 2a 2 uj 2 cos(0) - 2a 2 uj{t) 2 ) 

+ D(0) (-a 2 uj% sm{<p) - 2aw 3 ,w(i) cos(0) + 2auj 2 uj(t) 2 ) . 

Averaging over the fast oscillation, sin(</>) = leads to: 

*<«>-*<K i+ (£)')-^(£) ,+ ^^ 

The other extreme condition is the sudden limit when a > w. The fast variables 
can be integrated first leading to 

* N w-Ki-S i;S )(!)«» ™ 

and D staying constant up to first order in — . The dynamics can be translated to the 
final expectation at the end of the adiabat: 



H \ if i + m 2 i-(%y \ ( h 



l I - 2 I l-(^) 2 l + (^) 2 I \ L , 

' final v v <»i ' v «>i ' ' \ 7 initial 



(34) 



where u>i/f is either tOh or oj c depending on whether the evolution takes place on the 
power or on the compression adiabat. 

The dynamics of the working fluid on the different branches of the engine can be 

summarized by a branch propagator Uk which maps the set of operators A across the 

branch. An example is the hot isochore A(rh) — UhA(0). 
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Figure 1. The limit cycle of the engine in the (Se A 1 ) plane. The finite 
time cycle is indicated by the points A B C D. The cycle with infinite time 
allocation to all branches is indicated by the points W E Q F. Another long 
finite time cycle is indicated in yellow. The equilibrium isotherms of the hot 
bath intersecting at point E and the cold bath intersecting at point F, are also 
shown, fruih = 2, hu>c = 1, ksT^ = 5, ksT c = 1, T = 0.03 in units of u) c . 
Th = 6 , r c = 12 , "r>j C = = 1 in units of uj^ 1 - 



5. The Engine In Action 

We can now turn to look at the engine cycle as a whole. The four branches of the 
engine can be summarized as follows: 

(i) The hot isochore: A — ► B in figure ^ characterized by iOh the working medium 
frequency, the temperature of the hot bath, the heat transfer rate, 
(H) eq (Th) the equilibrium energy of the working fluid point E. The external 
control is the time allocation which determines the extent of equilibration. 
The dynamics of the operator set is represented by the propagator Uh ■ 

(ii) The power adiabat: B — ► C in figure ^ characterized by the change from u>h to 
u) c of the working medium frequency. The external control is the time allocation 

translated in to the adiabatic parameter cthc- The dynamics of the operator 
set is represented by the propagator Uh c - 

(iii) The cold isochore: C — > D in figure. ^ characterized by the uj c the working 
medium frequency, the temperature T c of the hot bath, r c the heat transfer 
rate, {H) eq (T c ) the equilibrium energy of the working fluid point F. The external 
control is the time allocation r c which determines the extent of equilibration. The 
dynamics of the operator set is represented by the propagator IA C . 

(iv) The compression adiabat: D — > A in figure ^ characterized by the change from 
oj c to ujh of the working medium frequency. The external control is the time 
allocation r c h translated to the adiabatic parameter a c h- The dynamics of the 
operator set is represented by the propagator U c h- 
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Figure 2. A limit cycle trajectory in the (H), (L) and (D) coordinates. A — > B 
is the hot isochore, B — » C is the power adiabat, C — » D is the cold isochore and 
D — ► A is the expansion adiabat. 



5.1. The limit cycle 

Once the engine is ignited it is expected that it will settle down to a smooth mode of 
operation determined by the control parameters. This asymptotic cycle will be termed 
the limit cycle. 

The dynamics of the four strokes of the cycle can be represented by the cycle 
propagator: 

Ucyc = Uch U c Uhc Uh, ■ (35) 

It is therefore expected that each point of the limit cycle trajectory is invariant to 

the propagator U C ycMimit = Mimit (we choose point A). The cycle propagator is a 
product of the branch propagators. They belong to the class of completely positive 
maps, U C yc therefore, is also a completely positive map. It has been proven that if 

Aumit is unique then a repeated application of U cyc will lead monotonically to the 
limit cycle [321 EH ■ This means that the limit cycle has to close. Consequently all 
variables are cyclic. This property allows to study the optimal performance of the 
operating engine. 

A typical limit cycle trajectory is shown in figure [21 The oscillation of (L) and 
(D) on the isochores with a different frequency Uh or w c is clearly seen. The canonical 
invariance condition is preserved by the limit cycle. Moreover even if the initial state 
of the system does not belong to the class of canonical invariant states after a few 
cycles it will approach monotonically such a state |1U| . 
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5. 2. Optimizing the operation of the engine 

Finite time thermodynamics allows us to optimize power as well as efficiency. We 
will assume, for simplicity, that the external parameters are fixed. These include 
the bath parameters Th and T c as well as the heat transport rates Th and r c . At 
first the frequencies u)h and lo c and with them the compression ratio C = u>h/u> c are 
also set. The optimization is carried out with respect to the time allocations on 
each of the engine's branches: Th, The, t c , and r c h- This sets the total cycle time 

T = T h + The +T C + T ch . 

The maximum efficiency limit is obtained for infinite time allocations on all 
branches. This cycle maximizes the work output. Its trajectory passes through the 
points W E Q F in figure ^ an d the efficiency reaches the value r) = 1 — Since 
for this limit the cycle time t is infinite, obviously the power of this cycle is zero. 

Optimizing the power is far more interesting, ft requires a finite cycle time 
t. As a result the optimizer has to compromise and shorten the time allocated to 
equilibration on the isochores. In addition restricting the time on the adiabats results 
in nonadiabatic conditions that expose the nature of friction. 

5.2.1. Optimizing the operation of the engine for quasistatic conditions. The power 
optimization in the quasistatic regime is determined by heat transport from the hot 
and cold baths. To reach quasistatic conditions, ample time is allocated to the adiabats 
(Cf. equation (J20J). As a consequence, the Number operator becomes invariant on 
the adiabats. 

The time allocations on the isochores determine the change in the number 
operator N = (N) (Cf. equation (J2SJ): iV B = e - ThTh (N A - N^ q ) + N^ q on the hot 
isochore, where iV B is the number expectation value at the end of the hot isochore, 
N A at the be ginning and AT™ is the equilibrium value point E. A similar expression 
exists for the cold isochore. 

The work in the limit cycle becomes 

W q = E C -E B + E A -E D = h(u c - uj h )(N B - iV D ) (36) 

where the convention of the sign of the work for a working engine is negative, in 
correspondence with Callen [S]. 

The heat transport from the hot bath becomes 

Q h = E B -E D = ^(N* - N D ) . (37) 

The efficiency rj becomes independent of the time allocations: 

r) q = Vo = 1 - — = 1 - \ ■ (38) 
uj h C 

In the limit cycle for adiabatic conditions N 13 = N A which leads to the relation 
n b = _Vf H—(N h - N c ) + N h (39) 

JV i _ e T c T c +T h T h V V eg ^eq) T ^ e q ■ 

By demanding a closed cycle the number change iV B — iV D is calculated, leading to: 
W q = h(u h - u c ){N* q - ^)L__iL__J (40) 

= - G w (T c ,uJc,T h ,u} h )F(x c ,x h ) , 

where the scaled time allocations are defined x c = r c r c , and Xh = T/jT^. The work 
W q lj^U|l becomes a product of two functions: Gw which is a function of the static 
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constraints of the engine and F which describes the heat transport on the isochores. 
Explicitly the function Gw is 

G w (T c ,u c ,T h ,u>h) = ^ (uj h - uj c ) ( coth ( J^t, ) - coth 



2 v ' V \2k B T h ) \2k B T c 

< F < 1, therefore for the engine to produce work Gw > 0. The first term in 
<|41[1 is positive. Therefore Gw > translates in to the restriction that ^ > ^ or 
in terms of the compression ratio 1 < C < ■p L . This is equivalent to the statement 
that the maximum efficiency of the Otto cycle is smaller than the Carnot efficiency 
Vo < Ve- 
in the high temperature limit when -jqj, <C 1, Gw simplifies to 

G w = k B T c {\ - C) + k B T h (l - C- 1 ) . (42) 
In this case the work W q = —Gw F can be optimized with respect to the compression 
ratio C — ^ for fixed bath temperatures. The optimum is found at C = \J~^- As a 
result the efficiency at maximum power for high temperatures becomes 

% = i - sjY h > ( 43 ) 

which is the well known efficiency at maximum power of an endoreversible engine 
El 123 ESI 133 • Note that these results indicate greater validity to the Curzon- 
Ahlborn result from what their original derivation % indicates. 

The function F defined in 14U|) characterizes the heat transport to the working 
medium. As expected F maximizes when infinite time is allocated to the isochores. 
The optimal partitioning of the time allocation between the hot and cold isochores is 
obtained when: 

r h (cosh(r c r c )-i) = r c (cosh(r h T ft )-i) . (44) 

If (and only if) = T c the optimal time allocations on the isochores becomes t% = t c . 

Optimizing the total cycle power output V is equivalent to optimizing F j t since 
Gw is determined by the engine's external constraints. The total time allocation 
t = Ti SO + T a di is partitioned to the time on the adiabats T a di which is limited by the 
adiabatic condition, and the time Ti SO allocated to the isochores. 

Optimizing the time allocation on the isochores subject to (|44|l leads to the 
optimal condition 

r c r(cosh(r, l r/ l )-l) = sinh(r, 1 T/ l +r c r c )-sinh(r c r c )-sinh(r, i r, l ) .(45) 

When Th — T c = V this expression simplifies to: 

2x + TT adl = 2 sinh(x) (46) 

(where x — T c t c = T/jT/J. For small x equation i ffifll can be solved leading to the 

i 

optimal time allocation on the isochores: t c — t/j w (^T ac ii/3) 3 /T. Taking into 
consideration the restriction on the adiabatic condition this time can be estimated 
to be: t c = Th ~ y (5) 3 ■ When the heat transport rate T is sufficiently large, the 
optimal power conditions lead to the Bang-Bang solution where vanishingly small time 
is allocated to all branches of the engine [62] . 

For the conditions r = T c = , and restricting the time allocation on the adiabat 
to be T a di = -}- + — ! the quasistatic optimal power production as a function of cycle 
time r becomes: 



V q = h(uj h - lu c ) (coth (gj^) - coth ( 



2k B T c 



tanh( ^ (r- 
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Figure 3. The cycle of the engine in the (Se ,T) plane (solid lines) and the 
(Svn iT) plane (dashed lines). The finite time cycle is indicated by the points 
A B C D. The difference between Se and Syjv decreases at the end of the 
isochores. The parameters used were: hui^ = 2, ksT^ = 5 hui c = 1, ksTc = 1, 
r = 0.05, Th = t c = 10 , t/ ic = Th c = 3. Notice that Se > Svn an d the non 
monotonic change in Se on the adiabats. 



which is an upper limit to the power production Cf . figures |5j |5] and [7| 

The entropy production ASu (|21|l reflects the irreversible character of the engine. 
In quasistatic conditions the irreversibility is completely associated with the heat 
transport. ASu can also be factorized to a product of two functions: 

AS U = G s (T c ,u! c> T h ,u> h )F(x c ,x h ) , (48) 

where F is identical to the F function defined in l|4U|) . The function Gs becomes: 

G s {T c ,u c ,T M ) = - ( - ( coth ( - coth ( Jf^L-)) .(49) 

K ' 2\k B T h k B Tj \ \2k B T c ) \2k B T h JJ V ; 

Due to the common F{x c ,Xh) function the entropy production has the same 
dependence on the time allocations Th and r c as the work W [S]. As a consequence 
maximizing the power will also maximize the entropy production rate ASu/t. Note 
that entropy production is always positive, even for cycles that produce no work as 
their compression ratio C is too large, which is a statement of the second law. 

The dependence of the G s function on the compression ratio can be simplified in 
the high temperature limit leading to: 

Gs = , (50) 

J- h ±c 

which is a monotonic decreasing function in the range 1 < C < ^jr- that reaches a 
minimum at the Carnot boundary when C = 

5.2.2. Lost work required to counter the friction Quasistatic conditions only obtain 
for infinitely long adiabats. To determine the validity and applicability of this limit, 
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we must consider shorter adiabats and verify our results to higher orders. We must 
consider finite adiabatic branches which will lead to friction (Cf. equation J2HJ0- 
The accumulated work required to overcome the friction is defined as: 

W f (t) = - f a(L)dt , (51) 
Jo 

for each adiabatic branch. To second order in the nonadiabatic coupling — , the 
accumulated work becomes: 

W f (t)= - a (J_<L)(O)Bin(0)- ^a( J ff)(Q) + i(D)(Q)) S in 2 (^) (f,2 
Averaging on the fast oscillating phase <\> leads to: 

2 

W/w£ff(H>(O) + |0>(O) (53) 

When (D)(0) = and T — > oo, the accumulated work to counter friction on both 
adiabats for the limit cycle becomes 

which at the high temperature limit simplifies to 




Equations (|54|) and l|55|l reflect the finite time allocation on the adiabat (i.e. the 
finiteness of a) which always reduces the produced work relative to the quasistatic 
case [5]. Remember that for an engine producing work W = W q + Wf should be 
negative. This justifies the use of the term "friction" . When the friction is sufficiently 
large, no work can be produced at all in the given engine. This can be seen in figure 
showing the quasistatic work W q corrected by Wf as a function of the temperature 
ratio r = 

The entropy production due to friction is 

As f = ; i [ - \ J^coth (J^A + f^V ^coth ( h - 




k B T h \2k B T c J \LUh J k B T c \2k B T h 

At the limit of high temperature (2J simplifies to AS f = ± (j^j 2 (c^ +C" 3 ^-V 

The nonadiabatic nature of the motion leads to an increase in entropy production in 
the baths, in accordance with our conceptions of "friction" . 

5.2.3. Optimizing the engine in the sudden limit The extreme case of the performance 
of an engine with zero time allocation on the adiabats is dominated by the frictional 
terms. These terms arise from the inability of the working medium to follow 
adiabatically the external change in potential. A closed form expression for the sudden 
limit can be derived based on the adiabatic branch propagator U^c and U c h l|34|l . The 
expressions are complex and therefore will not be presented. 

Insight on the sudden limit can be obtained when the heat conductance terms T 
become infinite. The performance therefore is solely determined by the friction. For 
this limiting case the work per cycle becomes: 

W s = 7iw c coth(— — - huj h coth — — ) .(57) 

4w c w/j V k B T h k B T c I 
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The maximum produced work — W s can be optimized with respect to the compression 
ratio C. At the high temperature limit: 



1. _ ,„ 2 ., t T r 1 



W s = ^k B T h {C z - l)(^f - j^) ■ (58) 



2 " " v "T h C 2 

For the quasistatic optimal compression ratio C — W s is zero. The optimal 

Y c 

compression ratio for the sudden limit is C = ( ^ ) , leading to the maximal work 



in the high temperature limit 




(59) 



The efficiency at the maximal work point is 

i - sfzi 

(60) 

Equation i|6U|) leads to the following hierarchy of the engines maximum work 
efficiencies: 

Vs < Vq < Vc ■ (61) 

Equation Ijtjlf) leads to the interpretation that when the engine is constrained by 
friction its efficiency is smaller than the endoreversible efficiency where the engine is 
constrained by heat transport that is smaller than the ideal Carnot efficiency. At the 
limit of T c — > 0, i] s = ~ and r\ q = r) c = 1. 

An upper limit to the work invested in friction Wf is obtained by subtracting 
the maximum work in the quasistatic limit 14U|) from the maximum work in the 
sudden limit l|57|l . In both these cases infinite heat conductance is assumed leading to 
N B — N^ q and N D — N^ q . Then the upper limit of work invested to counter friction 
becomes: 

W/ = Rw/, ^ ■ ( 62 ) 

At high temperature l|t)2|) simplifies to: 

1 k b T h (C-l) 2 (C- 2 + ^ 

* J-h 

Figure 01 compares the maximum produced work at the high temperature limit of the 
quasistatic and sudden limits for the compression ratio that optimizes the quasistatic 

work C = ■Jljp- and the sudden work C — (^i-) 1 / 4 . The work against friction W/ 

(equation (|63(l ) can be identified the figure as an increasing function of the temperature 
ratio. For the compression ratio that optimizes the quasistatic limit the sudden work 
is zero. This means that all the useful work is balanced by the work against friction 
Wf = W q . Zero produced work is obtained in the quasistatic limit for the compression 
ratio C = ^ and becomes negative for the sudden limit. 



W f = -k b T h (C - 1) 2 (C- 2 + ^) . (63) 
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Figure 4. Maximum produced work for the quasistatic — W 3 (magneta), 
quasistatic corrected by friction — W q — Wt, when — = 0.5 (purple) and 
sudden — W s (green) normalized by ksTh as a function of the temperature 
ratio r = ^jp-. The plot is at the high temperature limit for three compression 

ratios: C = r (triangles) corresponding to the Carnot limit, C = r 1 / 2 (solid) 
corresponding to the endoreversible optimum compression ratio and C = r 1 / 4 
(dotted) corresponding to the sudden optimum. The work against friction in the 



sudden limit W 



/ 



W — W s is indicated. 



5.3. Optimal Power Production 

The full range of the performance options of the engine is explored by finding the limit 
cycle numerically for a random set of time allocations Th, Thc> Tci T c hi on the different 
branches. 

Figure [S] displays the power as a function of total cycle time and efficiency 
for the compression ratio C = yjTh/T c that optimizes the quasistatic limit at high 
temperatures. All plots show a clear global maximum in efficiency at a specific cycle 
time. In all cases, the deviation of the maximum power point from the quasistatic 
estimate is not large. Therefore the efficiency at maximum power is quite close to 
the maximum Otto efficiency r\ Q At large cycle times the optimal power obtained 
by optimizing the time allocations fits very nicely the quasistatic limit. A closer 
examination of the cycles with optimal time allocation shows that they form a 
discontinuous set. At different total cycle times r different combinations of time 
allocation lead to maximum power. Underneath the optimal time allocation we can 
identify cycles which are classified as being sudden cither on the hot-to-cold adiabat 
or cold-to-hot adiabat. When the heat conductance T is high these sudden cycles 
outperform the quasistatic ones at short times. Cycles which could be classified as 
sudden on both adiabats are absent in this plot since they produce almost zero power. 
This is expected for the compression ratio C = \/Th/T c for which W s is zero (Cf. 
equation (ff)3|) ). 

Figure displays the power as a function of total cycle time and efficiency for 
intermediate temperature and a lower compression ratio of C = (^ L ) - 3 . Comparing 
cycles operating at high temperature (figure [SJ to ones at intermediate temperature 
(figure EJ) show significant similarities. The performance of the engine approaches the 
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Figure 5. The Power production as a function of the total cycle time (left 
panels) and as a function of efficiency (right panels) for the high temperature 
limit , ^rfi = 0.05 and fi^S = 0.1 (hui c = 0.025). The compression ratio is 

chosen to correspond to the optimal value in the quasistatic limit C = y^T^/Tc. 
Upper panels are for large values of heat transport (r = 0.6) and the lower panels 
are for small value of heat transport (T = 0.03). The graph was produced from 
a random sample of time allocations on all branches (gray filled squares). The 
black circles represent the optimal time allocation values which correlate to the 
quasistatic limit. The blue triangles classify the sudden limit on the cold-to-hot 
adiabat. The red squares classify the sudden limit on the hot-to-cold adiabat. 
The Carnot efficiency r) c , endoreversible (CA) efficiency r\ q the Otto efficiency r\ a 
and the sudden efficiency r) s are indicated in the right panels. The purple line 
represents V q , the quasistatic power limit 1471 . 



quasistatic limit when the heat conductance is small. A careful examination of figure 
show some cycles with low power output that exceed the Otto efficiency. This is 
not a contradiction since the limit restricts only quasistatic cycles. Thus some sudden 
like cycles can exceed this limit. In addition at high heat conductivity a reminiscent 
of some bang-bang cycles at very short cycle times can be identified. 

Operating the engine in the compression ratio corresponding to the optimum for 
the sudden limit C = (Th/Tc) 1 / 4 unravels a different set of optimal time allocations 
leading to maximum power. This can be observed in figure0which shows the power as 
a function of cycle time and its efficiency for this compression ratio. For conditions of 
high heat conductance the sudden cycles on both adiabats out-perform all other cycles. 
These cycles are very close to the bang-bang solutions. The efficiency at maximum 
power also corresponds to this limit (Cf. equation (||)UJ|). 
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Figure 6. The Power production as a function of the total cycle time (left panels) 
and as a function of efficiency (right panels) for the intermediate temperature limit 
jjTTj^ = 0-4 and = ^-0 (hui c = 1)- The compression ratio is chosen to be 

C = (^ L ) ' 3 . Upper panels are for large values of heat transport (T = 0.6) and 
the lower panels are for small value of heat transport (V = 0.03). The graph was 
produced from a random sample of time allocations on all branches (gray filled 
squares). The black circles represent the optimal time allocation values which 
correlate to the quasistatic limit. The blue triangles classify time allocations 
with a shorter time on the isochore relative to the black circle case. The Carnot 
efficiency r] c , endoreversible (CA) efficiency r) q the Otto efficiency r) and the 
sudden efficiency rj s are indicated in the right panels. The purple line represents 
the quasistatic power limit 1471 . 



6. Discussion 

The performance of the model Otto engine shows close resemblance to macroscopic 
real engines. The power has a definite maximum which is the result of a compromise 
between minimizing heat transport losses and frictional losses |35j . How does this 
phenomena emerge from the quantum origins of the engine? Put more boldly what is 
quantum in the engine? 

The quantum laws governing the heat transport from the baths to the working 
medium (equations (j3J| and (|22|l ) have been derived from first principles (231221 m 
the Markovian limit. This limit guarantees the tensor product form of the extended 
system, allowing for a separation between the local variables of the system and the rest. 
This separation is a prerequisite for a thermodynamic description |3b| . Nonmarkovian 
models may also lead to a similar structure [23, but are beyond the scope of the 
present study. The irreversibility in these equations can be interpreted as emerging 
from a quantum measurement 1201 1381 . Specifically for the present system they 
coincide with the linear Newtonian heat transport law at high temperature. 
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Figure 7. The Power production as a function of the total cycle time (left 
panels) and as a function of efficiency (right panels) for the high temperature 
limit jjpy - = 0.02 and ^ = 0.16 (hw c = 0.1). The compression ratio is chosen 

to correspond to the optimal value in the sudden limit C = ^ • Upper 

panels are for large values of heat transport (r = 0.6) and the lower panels are 
for small value of heat transport (T = 0.03). The graph was produced from 
a random sample of time allocations on all branches (gray filled squares). The 
black circles represent the optimal time allocation values which correlate to the 
quasistatic limit. The blue triangles classify the sudden limit on the cold-to-hot 
adiabat. The red squares classify the sudden limit on the hot-to-cold adiabat. 
The Carnot efficiency r) c , endoreversible (CA) efficiency r\ q the Otto efficiency r\ 
and the sudden efficiency r) s are indicated in the right panels. The purple line 
represents the quasistatic power limit 1471 . 



The quantum nature of the heat transport on the isochores becomes significant 
at low temperatures when > 1. Planck's constant relates the frequency to a 
discrete energy unit Huj. The expressions for work (equations (|4*U1) and J57|)) and 
power production (equation (|4"7ll ) are linearly related to this Quant of energy. At 
higher temperatures the energy unit is replaced by fcgT and h disappears. Moreover 
the thermodynamical expressions for maximum power and efficiency coincide with well 
established relations known from finite time thermodynamics jjH HQ] . The finding is 
in line with the study of Geva et. al. [7j which has shown that although the quantum 
heat transport laws are not linear in general (equation l|24|) '). and the system can be far 
from equilibrium, nevertheless they lead at high temperatures to the Curzon Ahlborn 
efficiency at maximum power. 

The rapid dynamics on the adiabatic branches is responsible for the frictional 
behaviour. The quantum origin can be traced to the noncommutability of the 
Hamiltonian at different times. As a result the system is not able to adiabatically 
follow the instantaneous energy representation. For the harmonic model the origin of 



Irreversible performance of a quantum harmonic heat engine 



22 



the noncommutability comes about from the kinetic energy not commuting with the 
potential energy, the origin of which is the fundamental relation [Q, P] = ih. Plank's 
constant appearing in the commutation relation of the kinetic and potential energy 

^%-\ = ih~H}-D (Cf. equation JJJ) is cancelled by the inverse of Planck's constant 
appearing in the equation of motion or Ij22(l . 

The deviation from adiabaticity is independent of the mass of the oscillator. This 
inertial effect appears only in the combination a — uj/lo the nonadiabatic parameter. 
This means that friction like behaviour is expected whenever the internal timescale of 
the working medium is comparable to the engines cycle time. Without a dissipative 
branch, however, this "friction" cannot be interpreted as such - it is even reversible. 

The most important noncommutability of the engine is between the branch 
propagators %i c hMc\ 7^ or ^AhcMh] 7^ 0; composing the cycle propagator U cyc in 
equation (|35|l . The sequence of alternating propagations of adiabats and isochores 
is what breaks the time reversal symmetry and allows the engine to produce power 
irreversibly. 

7. Conclusion 

The quantum Otto cycle model engine served its purpose of unravelling the emergence 
of thermodynamics out of quantum mechanics. Due to the simplicity of a working 
medium composed of harmonic oscillators the following insights have been illustrated: 

• Thermodynamical variables are related to quantum observables which are 
sufficient to reconstruct the state of the working medium. The energy variable 
is supplemented by additional variables required to describe the deviation from 
adiabatic behaviour in the limit cycle. The knowledge of these observables enables 
the calculation of thermodynamical variables such as entropy and temperature. 
The first law emerges as a statement accounting for the energy change in an open 
quantum system. The second law is related to the fact that the energy entropy is 
always larger than the von Neumann entropy which is invariant to free evolution. 
The information required to restore the von Neumann entropy from the energy 
entropy is partially erased on the isochores. The open dynamics on the isochores 
can be considered as a partial measurement of energy. As a result information on 
the operators L and D which do not commute with energy is lost. 

• Quantum origins of irreversible losses are due to heat transport and with it the 
unavoidable dephasing. Quasistatic cycles which almost follow adiabatically the 
energy frame show irreversible losses dominated by heat transport. The optimal 
performance of these cycles closely resembles the endoreversible limit. 

• Deviation from quantum adiabatic behaviour is reflected by frictional losses which 
emerge from the generation of inertial components on the adiabats and their 
dephasing on the isochores. The frictional loss terms are independent of the cycle 
direction of motion. 

• At high temperatures, explicit quantum variables and constants disappear. The 
performance of the engine shows universal behaviour which is a function of 
temperature T compression ratio C and the ratio between the cycle frequency 
proportional to a and the internal frequency of the working medium w. 
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Appendix A. Canonical invariance for the product form 

The canonical invariance of a density operator represented as a product of exponents 
is considered: 

p = e^X 1 e^ 2 * 2 ...e> Mhjtk ...e> MNjtN (A.l) 

where the set {X} form a closed Lie algebra and p = p(t) are time dependent 
coefficients. Canonical invariance means that the form of l|A.l|) is preserved during the 

evolution. To test this hypothesis it can be noted that pp -1 forms a linear combination 
of the set of operators {X}: 

^-p- 1 = p 1 X 1 + p 2 e^X 2 e-^ + pze^e^X^^e-^* 1 + ...(A.2) 
at 

The Baker- Hausdorff relation e A Be~ A = B + [A,B] + \[A, [A, B]]... and the 
closure of the Lie algebra to commutation means that the RHS of equation (|A.2J| is a 
linear combination of the operators {X} in the algebra. 

For canonical invariance to prevail during the evolution also should 

become a linear combination of the set of operators {X}. This will lead to N coupled 
differential equations for the coefficients p(t). 

The simpler case is when the Hamiltonian is also composed of a linear combination 
of the same algebra then = — ih~[H, p]/5 _1 = —ifoH — ihpHp^ 1 forming a linear 

sum of the same set {X } and canonical invariance prevails. 

For an open quantum system the equivalent statement would be that in addition 
Coip)^ 1 also forms a linear combination of the set {X}. This condition becomes: 

CoipYp- 1 = £ (^SPP}r 1 -\(t'}F j +pF j F^- 1 )j (A.3) 

where F are Lindblad operators [19; . Equation l|A.3|l will form a linear combination 
of the set of operators {X} if the operators FjPj and FjpPJ p~ x are part of the set 
{X}. 

Specifically for the harmonic oscillator set P 2 , Q 2 and D = (QP + PQ) and the 
Lindblad operator F — a this condition is met. On the contrary for F — H , H 2 is 
not part of the Lie algebra therefore, canonical invariance is not met in the case of 
pure dephasing. 
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Appendix B. The relations between the parameters in the product form 
and the exponential form 

The state of the system can be expressed in an exponent product form HlOjl or as 
an exponent of an operator sum 117|) . To determine this equivalence, the expression 
pap -1 and pa) p^ 1 are evaluated. The normalization Z is identical for both forms 
and therefore is not considered. In the exponent-product form IjlOl) using the Baiker- 
HousdorfF formula one obtains: 



pap 



-1 _ (Bhw 



( e /3^ _ Ae~^-fj*)a - 2e- /3Sw 7*at (B.l) 



ptfp- 1 = 2e- 0hw ja + e-^at (B.2) 
Using the operator sum form (|17|) to evaluate the same two expressions, leads to: 



pap 1 



cosh(^/ {hujX2f - 4xixl) - X2 , sinh(y / (hw X 2) 2 ~ 4Xix!) ) a 

\l (hcJX2) 2 - 4xiX* 

2Xj = sinh(vW2) 2 -4xiXt) I ^ (B.3) 



'(hu;x2) 2 -4xiX* 



and 



Tied 



patp 1 = [ cosh(y (Hujx2) 2 - ^XiXi) + X2 , = sinh(y {fiuj X 2) 2 - 4xiXi) ) a f 

V ( huj X2) 2 - 4xiXi 



2xi 



,\J{T^X2) 2 -4xiX* 
Equating both formula pairs results in the identification 



sinh( v /(^x 2 ) 2 -4xiXi) a (B.4) 



arcsinh Uje"^ ^I r+^-i-Urf 



Xi = 2 v ; '- (B.5) 

'-1677* + (e 2 ' 3 ' i ^ ) 2 



7" 



7* 1 - e 2/3ri " + 477* 



Xi = Xi— , X2 = Xi" % 
7 Zjnuj 

which establishes the relation between 7 and f3 in the product form (|17|) and x in the 
the operator sum form i|l(J|) . 
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